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Abstract. Filaments on the surface of a microorganism such as Paramecium or Ophalina beat highly syn- 
chronized and form so-called metachronal waves that travel along the surfaces. In order to study under 
what principal conditions these waves form, we introduce a chain of beads, called rowers, each periodically 
driven by an external force on a straight line segment. To implement hydrodynamic interactions between 
the beads, they are considered point-like. Two beads synchronize in antiphase or in phase depending on 
the positive or negative curvature of their driving-force potential. Concentrating on in-phase synchroniz- 
ing rowers, we find that they display only transient synchronization in a bulk fluid. On the other hand, 
metachronal waves with wavelengths of 7-10 rower distances emerge, when we restrict the range of hydro- 
dynamic interactions either artificially to nearest neighbors or by the presence of a bounding surface as in 
any relevant biological system. 

PACS. PACS-key discribing text of that key - PACS-key discribing text of that key 



1 Introduction 

For microorganisms with sizes between 1 to 1000/^m the 
surrounding fluid appears to be very viscous. They, there- 
fore, need a highly specialized locomotion apparatus in 
order to propel themselves [T1[2J|21I1]. Paramecium, e.g., 
uses a dense array of whip-like appendages on its surface, 
called cilia, that move fluid close to its surface by beat- 
ing periodically with a succession of power and recovery 
strokes [S]. The properties of cilia have been studied for 
decades [6]. In particular, the coordinated beat cycles of a 
multitude of them, observable as traveling waves [7|, has 
received a lot of attention recently IHUllIinillllinillllllllSl 
[16]. The formation of traveling or so-called metachronal 
waves needs a high degree of synchronization between the 
beating cilia. Many different types of these waves called 
sym-, anti-, laeo- or dexioplectic arc found on a variety 
of ciliated microorganisms, on certain types of comb jelly, 
and in the respiratory tract of mammals [T7, 18, 19, 20, 21, 

Synchronization in large populations of interacting dy- 
namic elements can be observed in many physical, chem- 
ical, biological, as well as social systems and a unifying 
description was formulated by Kuramoto f23| . Examples 
range from micromechanic resonators used for the con- 
struction of highly sensitive mass-, spin-, and charge- mea- 
suring devices P31|251I26] , over the self-organized beat cy- 
cle of single cilia, where the action of dynein motors is 
coordinated by the curvature of the cilium |27j , the syn- 
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chronized flashing of fireflies, the working of muscles and 
neurons, to the synchronization of the applause by an au- 
dience after the last bar of a classic concert [281129] . 

In particular, in viscous fluids hydrodynamic coupling 
leads to the synchronization of rotating helical filaments 
such as bacterial flagella [501I3T] , rotating paddles [32] , mi- 
crofluidic rotors [551IM1I55] , the pair of beating flagella in 
Chlamydomonas |36j , or flagella of neighboring sperm cells 
[571155113511^ . The most astonishing example for synchro- 
nization in viscous fluids are metachronal waves. A lot of 
works have been devoted towards understanding their for- 
mation. They either use minimal models that abstract cilia 
as point-like objects in order to identify essential features 
[5M [TUl[m[T^ or are based on a more accurate modeling of 
the beating cilia [T51[Hl[T51ll6| . From all these studies the 
common picture has evolved that hydrodynamic interac- 
tions can synchronize beating cilia and thereby cause the 
formation of metachronal waves. 

This work uses the minimal model of Coscntino Lago- 
marsino et al. [5] that is based on a linear array of driven 
oscillators, so-called rowers, to clarify under what condi- 
tions metachronal waves are able to form. First, we will 
demonstrate that the long-range nature of hydrodynamic 
interactions impede the formation of metachronal waves. 
Only when they become short-ranged close to bounding 
surfaces do metachronal waves occur. Secondly, a stroke 
that becomes faster during one half cycle leads to in-phase 
synchronization of a pair of oscillators and thereby to 
metachronal waves in a linear array of rowers with rel- 
atively long wavelengths. Clever recent experiments of a 
pair of driven colloids have demonstrated antiphase syn- 
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chronization when the stroke slows down during one half 
cycle [12]. It is our hope that these experiments can also 
demonstrate in-phase synchronization. 

The paper is organized as follows. In Sec. [2] we intro- 
duce our model including the potential for the force to 
drive the rowers. In Sec. [3] we investigate the synchro- 
nization of two rowers by distinguishing between nega- 
tive and positive curvature in the driving-force potential. 
Then we study synchronization in a chain of rowers in 
Sec. m We first introduce an order parameter for iden- 
tifying metachronal waves, demonstrate their formation 
in a bulk fluid when we artificially restrict hydrodynamic 
interactions to nearest-neighbor rowers, and show that a 
bounding surface has the same effect. We finish with a 
conclusion. 




Fig. 1: Geometry of the open linear chain of rowers that 
move back and forth on line segments (dashed lines) with 
length 2s. The line segments are oriented along the y axis, 
their centers are separated by the lattice vector c. and the 
tilt angle between the segments and c is /3. 



2 The model 

In order to study basic features of metachronal waves, we 
abstract cilia as point-like beads following the successful 
works of Cosentino Lagomarsino et al. [5] , Niedermayer et 
al. [TT] , and Vilfan and Jiilicher [5] . We place the beads on 
a linear, periodic array and constrain their motions onto 
line segments, on which they move back and forth. This 
model was first introduced by Cosentino Lagomarsino et 
al. and we adopt their convention where a bead is called 
rower [8]. 

At low Reynolds number, the dynamics of the rowers 
is completely overdamped and the equations of motion for 
N driven beads read 

iV 

— ^ ^ rani,'^ van) n • (1) 
n=l 

Here is the position vector of bead m and Fm the force 
acting on it. The mobility matrices Mmn depend on the 
difference vector rmn = fm — fn- In an unbounded fluid 
with viscosity 77 and in the approximation of point-like 
particles, one uses 

{i^o^ ^ e-^ a if ri = m 

Sttti r \ r-^ / ^ ' 

where jio = 1/ {dnrja) is the Stokes mobility of a particle 
with radius a, G^^{r) the Oseen tensor, (g) means dyadic 
product, and r = \r\. In the course of this paper, we will 
also study the chain of rowers close to an infinitely ex- 
tended plane wall that bounds the fluid as depicted in Fig. 
[71 Then, the mobilities of point particles are more compli- 
cated. We summarize them in Appendix 1X1 In particular, 
the cross mobilities are given by the Blake tensor G™{r). 
In the following, we employ a rescaled version of the dy- 
namic equations ([T]) by introducing reduced quantities. 
We write the mobilities in units of /xq, the applied forces 
in units of a characteristic force value -Fq, and lengths in 
units of s, where 2s is the length of a rower's line segment. 
Then time is measured in units of s/IplqFq). 



Figure \T\ depicts the geometry of the linear chain of 
rowers separated by the lattice vector c. The chain has 
open ends for which we present most of our results. The 
effect of periodic boundary conditions is briefly discussed 
in Sec. 14.21 While Cosentino Lagomarsino et al. and Kotar 
et al. let the rowers beat parallel to the lattice vector c [S] 
[T2] , we introduce the tilt angle (3 of line segments against 
the lattice vector c as global parameter. The line segments 
arc oriented along the y axis and their centers are located 
at mc = mc(sin/3, cos /3). Together with the asymmetry 
of the power and recovery stroke introduced below, this 
tilt breaks the left-right symmetry of the linear chain of 
rowers and metachronal waves can then propagate with a 
definite direction either to the left or to the right. 

The driving forces acting on the rowers are always par- 
allel to the line segments. Furthermore, we only consider 
the y components of the hydrodynamic-interaction terms 
in Eq. ^ and neglect their x components. This reduces 
the spatial degrees of freedom of each rower to one and 
its position is given by = "mc + y-m^y (Fig. [1]). In or- 
der to check this approximation, we allowed excursions of 
the rowers perpendicular to their line segments by intro- 
ducing strong harmonic potentials along the x direction. 
In simulations of this two-dimensional rower model, we 
only found minor quantitative but no qualitative differ- 
ences to the results presented below in various examples. 
We therefore decided to work with the one-dimensional 
model, which needs much less simulation time. 

In addition to the continous displacement variable y„i, 
we use the discrete "geometric-switch" variable cr,„ = ±1 
to describe the rower's state. It was first introduced by 
Gueron and Levit-Gurevich [T3]. It determines the direc- 
tion of the driving force and hence if the rower moves 
parallel to Sy ((t,„ = 1) or antiparallcl to (dm = —1). 
The geometric switch changes sign when the rower reaches 
the displacement s from the center of the line segment. 
Therefore, the rower performs a periodic motion and can 
be regarded as a driven oscillator. For a driving force with 
constant magnitude Fq and non-interacting rowers, the 
period of this oscillation is r = As/ {iiqFq). In several fig- 
ures presented below we will refer time t to the period r. 
For later use we assign to each rower a phase variable that 
grows by 27r during one cycle starting from the center of 
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Fig. 2: Potential of the driving force for e = and e = —0.7 
and a = 0.2. The solid lines belong to cr = 1 and the 
dashed lines to cr = —1. For positive a the fast power 
stroke is performed from left to right. 
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Fig. 3: The synchronization of two rowers depends on the 
curvature e of the driving-force potential. The time evolu- 
tion of the relative phase Aip = (p2 — is shown for two 
intitial values, Aip = tt/2 and — 37r/4. Power and recovery 
stroke are symmetric, i.e. a = 0, the tilt angle /? = 7r/4, 
and the distance of the rowers is c = 0.28. 



the line segment = 0): 

fo if 2/™ G [ 0, 1) A (t™= 1 

TT I 

(fm = 27rn,„-|--cr„,2/„, + < tt if y,n G [ 1, -1) A <t„i = -1 
[27r if 2/„, G [-1, 0) A a„,= 1 

(3) 

Here, G N is increased by one after each completion of 
a cycle. 

We define the reduced driving force on the 771-th rower, 
Fmiym,crm) = ~dV {ym, cr^) / Oym, using the reduced po- 
tential 

V{y,n, CTm) = (^-(TrnVm + (acr„ + 1) + Vp . (4) 

The linear term on the right-hand side gives a constant 
force with magnitude Fq and the second term contributes 
a harmonic part to the potential with curvature e G (— 1, 1). 
The asymmetry parameter a G (—1,1) distinguishes be- 
tween a fast power and a slow recovery stroke. Curvature 
e and asymmetry a are the two important parameters of 
the potential. Finally, we also introduce a harmonic po- 
tential Vp centered at the switching points (y,„ = ±1). 
It acts only when the rower passes the switching point 
due to hydrodynamic drag forces from neighboring rowers 
and, thereby, prevents the rower to trespass the switch- 
ing point significantly. In Fig. [2] we illustrate the reduced 
potential for a = 0.2 and for two curvature values, e = 
and —0.7. The geometric switch variable is represented 
by solid (cr ~ 1) and dashed (cr = —1) lines, respectively. 
With the choice of a positive a = 0.2, the fast power stroke 
is performed at cr = 1 when the rower moves in positive y 
direction. 

Note that Cosentino Lagomarsino et al. [5] define the 
power and recovery stroke by introducing different Stokes 
mobilities in the forward (cr = 1) and backward (cr = 
— 1) motion. Such a choice does not directly affect the 
hydrodynamic-interaction terms in Eq. ([l} since the driv- 
ing forces are the same. We tested that this method docs 



not effectively break the left-right symmetry of the rower 
chain in order to create metachronal waves which run in 
one direction along the whole chain. In contrast, by al- 
tering the strength of the driving force we are able to 
break the directional symmetry of the chain and induce 
metachronal waves with a definite direction. 

A single rower performs now the following oscillatory 
motion in a potential with negative curvature as illus- 
trated in Fig. [31 It starts off slowly from a switching point 
(e.g. at ym = —1) and increases its velocity until it reaches 
the opposing switching point (e.g. at ym = 1). Here the 
switch variable a changes sign and the rower is dragged 
back into the opposite direction. Both, direction and mag- 
nitude of the velocity change discontinuously at the switch- 
ing points. The non-zero curvature e of the potential in- 
creases the period of a single rower relative to r, the period 
when only a constant force Fq acts. 



3 Synchronization of two rowers 

Figure [3] summarizes our numerical studies on the dynam- 
ics of two rowers. The time evolution of the relative phase 
Aip = ip2 — <Pi crucially depends on the curvature e of the 
driving-force potential. While positive e leads to antiphase 
synchronization regardless of the initial phase difference, 
the rowers assume the same phase when e is negative. For 
a constant force (e = 0), Aip keeps its initial value and 
does not change in time. We now develop an understand- 
ing for this behavior. 

For zero curvature e of potential V and identical power 
and recovery stroke, a = 0, the driving force on both par- 
ticles is constant and identical throughout one beating 
cycle. Then the mutual hydrodynamic interactions are al- 
ways equal in absolute value and thus the phase difference 
between both rowers does not change. The latter is even 
valid for the phase difference after one beating cycle if 
a 0. Due to the different forces in the power and recov- 
ery stroke, particles change their relative phase when they 
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move in opposite direction. However, during one beating 
cycle each particle performs both the power and the re- 
covery stroke. So, after one period the phase difference 
assumes its initial value and synchronization does not oc- 
cur. We stress again that a is significant only for breaking 
the left-right symmetry in a chain of rowers in conjunc- 
tion with the tilt angle /3 which then enables metachronal 
waves that travel in one direction only. 

A non-zero curvature in the driving-force potential breaks 
the symmetry in the hydrodynamic interactions between 
two rowers and leads to synchronization as already demon- 
strated in Fig. [21 Whilst Cosentino Lagomarsino et al. 
restrict e to be positive, which leads to an antiphase syn- 
chronization [51ll2j. we also allow negative values, thereby 
obtaining the desired in-phase synchronization. The sig- 
nificant part of the synchronization is accomplished when 
both rowers move in opposite directions. We first consider 
the case e < and assume that the preceding rower has 
just reached its switching point. After activating the geo- 
metrical switch, the preceding rower experiences a weaker 
driving force than the succeeding rower before the switch- 
ing point (see Fig. [5]) . Therefore, the flow field initiated by 
the succeeding rower is stronger and slows down the pre- 
ceding rower more strongly than vice versa. This leads to 
a decrease of the phase difference Aip until the succeeding 
rower also reaches the switching point. Then both row- 
ers move in the same direction and the phase difference 
Alp increases again. However, in total a decrease is real- 
ized during one cycle as illustrated in Fig. |21 For positive 
e, the preceding rower slows down the succeeding rower 
more strongly and Atp increases on average during one 
cycle until it reaches the value tt. To conclude, in the time 
evolution of the phase difference Alp = is a stable and 
Alp = ±7r is an unstable fixpoint for e < 0, whereas for 
positive e the stability of the fixpoints is reversed. 

In a very attractive experimental realization of the 
two-particle system, Kotar et al. use optical traps to cre- 
ate the driving-force potential with positive curvature e > 
[12]. Each colloidal bead is confined in the harmonic 
regime of the potential well created by a laser beam. The 
potential well switches between two positions such that 
a bead is always driven towards the potential minimum. 
The "geometric switch" is then initiated before the bead 
reaches the minimum. In agreement with our results, the 
two beads synchronize in antiphase. The trap potential 
created by an optical tweezer has a regime further away 
from the minimum, where the curvature is negative [42] . 
We propose to use this regime to realize in-phase synchro- 
nization experimentally. 

4 Synchronization in a chain of rowers 
4.1 Complex order parameter 

The collective dynamics of a chain of several hundred row- 
ers strongly depends on the parameters introduced above 
that determine the driving-force potential and the geom- 
etry of the rower chain. In order to classify the collec- 
tive dynamics and, in particular, to identify metachronal 



waves, we introduce a complex order parameter and study 
its longterm behavior. 

For an array of N oscillators, we can introduce N — 
1 phase differences . . . , Z\(/3iv-i} and their phasors 

z„ = e'^"^" G C, where Z\(p„ = </Jn+i — Vn- We define our 
complex order parameter as the mean of the phasors, 

N~l N-l 

N -I ^ N-l ^ ^ ' 

n—l n—1 

The magnitude A and polar angle 'P lie in the respective 
ranges A G [0, 1] and <P G [— tt, tt). We add a note here. In 
order to determine if oscillators fully synchronize to equal 
phases, usually the order paramter averages over the pha- 
sors e"^" for phases (pn P51l29j . However, this definition is 
not appropriate for identifying metachronal waves. 

For a large number TV of rowers with random phases, 
the phasors z„ are uniformly distributed over the unit cir- 
cle in the complex plane. Their mean, the complex order 
parameter Z, is close to 0, i.e., A sa with arbitrary 
^. If, on the other hand, the system displays a stable 
metachronism, where pairs of neighboring oscillators are 
phase locked with the same phase difference [55], then 
Zn = Zm for all phasors m,n G {1,...,N — 1}. In this 
case, Z is stationary. It lies on the unit circle in the com- 
plex plane = 1) and is equal to the phase difference 
of each pair of oscillators. 

In the following, we will also observe states where the 
long-time limit of Z fluctuates around a constant value. 
We therefore introduce the average values 

A=^ r A{t)dt ^=^,f mdt , (6) 

where T is the final time of the simulation run and At 
an appropriately chosen time interval. In addition, we also 
quantify fluctuations of the order parameter by calculating 
the standard deviations <ja and cr<j of the time averaged 
values. 



4.2 Synchronization in an unbounded fluid 

In Figs.|4|a) andSJ^b) we illustrate the temporal evolution 
of the order parameter Z{t) for an open chain of 200 rowers 
in different situations. For e = and a = 0.2 [golden graph 
(—-)], magnitude A{t) and, in particular, ^(t) strongly 
fluctuate in time. Their respective mean values are ^ ^ 1 
and <P ^ TT. Apart from a tendency of neighboring rowers 
to beat in antiphase, a noticable structure formation does 
not occur. 

We now introduce a negative curvature e = —0.3 in 
the driving-force potential. The green graphs (^) in Figs. 
lU^a) and SJ^b) still show a strongly fluctuating order pa- 
rameter. In contrast to e = though, we observe a sig- 
nificant increase of the time-averaged magnitude A. To- 
gether with <P ^ this indicates that some of the rowers 
synchronize in phase. A further understanding is achieved 
from Fig. [SJa) , where we plot the color-coded phases over 
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Fig. 4: Absolute value (a) and polar angle (b) of the order 
parameter Z{t) as a function of time for a chain of 200 
rowers for different parameter sets. The global parameters 
are c = 0.28 and /? = it/A and for all graphs the rowers 
start with the same random distribution of phases. The 
vertical red lines (^) mark the time inter vail At = 100 
for which we calculate the time averages A and ^. They are 
shown as dots on the right-hand side of the corresponding 
diagrams together with crrorbars indicating the respective 
standard deviations <ta and a^. Golden line (■— ): a = 0.2 
, e = 0; and long-range hydrodynamic interactions (l.r. 
HI) are used. Note for this parameter set we have added 
27r to all negative values of <P{t) to better represent its 
time dependence. Green line (— ): a = 0.2 and e = —0.3 
with l.r. HI. Magenta line ( — ): a = 0.2 and e = -0.3 
but hydrodynamic interactions are artificially restricted 
to nearest neighbors (n.n. HI). A metachronal wave occurs 
that travels in a definite direction. Blue line (— ): same 
situation but a = —0.2. The mean phase ^ is reversed and 
the metachronal wave travels in the opposite direction. 



time for the first 50 rowers. At t = we start with curva- 
ture e = and switch to e = —0.3 at t = 5. Then chain 
segments of approximately 8 rowers appear that are fully 
synchronized in phase as expected for negative curvature. 
However, after some time these clusters break apart and 
form somewhere else. The synchronization is only tran- 
sient since the rowers in the center of a cluster have re- 
duced friction coefficients compared to rowers at the edges 
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Fig. 5: Color coded phases of the the first 50 rowers in 
a chain of 200 rowers as a function of time t. The global 
parameters are a = 0.2, c — 0.28, and /3 = it/A. The 
simulations start with e = and at < = 5 the curvature 
changes to e — —0.3 in order to visualize the onset of 
synchronization, (a) For long-range hydrodynamic inter- 
actions transient in-phase synchronization occurs locally, 
(b) For nearest-neighbor interactions metachronal waves 
are visible. 



due to long-range hydrodynamic interactions. Hence, row- 
ers in the center move faster and destabilize the cluster. 

The reported instability is caused by the long-range 
nature of hydrodynamic interactions. Following Ref. 
we artificially make them short-ranged by including only 
the nearest neighbors in the sum of the dynamic equa- 
tions (H)). The range of hydrodynamic interactions can be 
controlled close to a surface, as we demonstrate in the 
next section. Now metachronal waves appear that travel 
from one end of the chain to the other with a wavelength 
of approximately eight rowers. Fig. [5jb) demonstrates im- 
pressively how the initially uncorrelated rowers evolve into 
a metachronal wave due to phase-locking. The process is 
completed for the whole chain after t = 175. This is also 
visible in Fig. Hl^a) [magenta line (— )] when A{t) « 1 is 
achieved. The phase difference [Fig.|4|b)] converges to ap- 
proximately «P ss 7r/4 which gives the corresponding wave- 
length of 8 rower distances. The very small standard devia- 
tions of A and ^ are a further measure that a metachronal 
wave has been formed nearly perfectly. Finally, reversing 
the sign of a [blue curves (— ) in Figs. Hlja) andlH^b)], 
which means exchanging the fast power and slow recovery 
strokes, also reverses the sign of 3 and the metachronal 
wave travels in the opposite direction. 

To break the left-right symmetry of the rower chain, 
parameter values a 7^ and (3 ^ -k/2 are needed but also 
the distance c of the rowers has to be sufficiently small 
for creating metachronal waves travelling in one direction. 
For a — or /3 = 7r/2 or by simply choosing c > 1, we 
observe metachronal waves running from the center of the 
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Fig. 6: Averaged magnitude A and polar angle 3 (inset) 
of the order parameter as a function of the number N of 
rowers in open (circles) and closed (squares) chains. 



chain into both directions. This is reminiscent to studies 
of Niedermayer et al. [TT] , where beads move on circular 
trajectories. 

Positive curvature e > of the driving-force poten- 
tial reproduces the results of Cosentino Lagomarsino et 
al. [5] . The tendency of neighboring rowers to synchronize 
in antiphase leads to a metachronism with shorter wave- 
lengths of approximatly 4 rowers. Compared to the case 
with e = —0.3, the system needs about ten times longer 
to reach the stationary state. 

We shortly discuss how the chain length of open chains 
and also periodic boundary conditions influence the oc- 
curence of metachronal waves. We concentrate on the same 
parameters as before (e = —0.3, /3 = 7r/4, and a = ±0.2) 
and use hydrodynamic interactions restricted to nearest 
neighbors. This allows us to realize the periodic bound- 
ary conditions by a closed chain where the first and last 
rower interact with each other. Figure |5] compares A and 
3 (inset) for open (circles) and closed (squares) chains for 
different chain lengths. For chain lengths above 50 rowers 
open and closed chains develop metachronal waves with 
nearly the same A and 3. So finite size effects become 
negligible. At shorter chain lengths the dynamics of open 
and closed chains differ strongly. Whereas for TV = 10 
(first data point) the same magnitude A ss 0.6 indicates 
some order but certainly no metachronism, the polar an- 
gles reveal different behavior in the two cases. For open 
chains ^ = 0, which means that central rowers synchro- 
nize in phase whereas the outer rowers lag behind due to 
the larger friction they experience. On the other hand, in 
closed chains, where all rowers have an identical environ- 
ment, a non-zero mean polar angle ^ is observed. A closer 
inspection of the individual rower phases shows that pairs 
of rowers almost fully synchronize in phase whereas they 
leave a larger phase gap to the next pair. We discuss these 
points so detailed here because in an experimental realiza- 
tion finite size effects of the chain may become important. 




Fig. 7: A chain of rowers a distance h above an infinitely 
extended wall. 
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Fig. 8: Inset: The relevant component of the Blake 
tensor as a function of h for, both, nearest-neighbor 
(|m — n| = 1) and next-nearest neighbor (|m — n| = 2) row- 
ers. The distance of the nearest-neighbor and next-nearest 
neighbor rowers are c = 0.28 and 2c = 0.56, respectively. 
The dashed lines indicate the respective values of the Os- 
een tensor, i.e., for height h — >■ oo. Large graph: ratio of 



both components, Gfjd"! — n\ — 1)/G22(| 



2). 



4.3 Synchronization near an infinitely extended wall 

So far we artificially restricted hydrodynamic interactions 
to nearest neighbors in order to make them short-ranged. 
In reality, the range of hydrodynamic interactions can be 
tuned by placing the chain of rowers a distance h above an 
infinitely extended wall. As sketched in Fig. [71 the rowers 
beat parallel to this wall. In order to describe the cross 
mobilities of point particles close to a wall, the Oseen ten- 
sor has to be replaced by the Blake tensor [JT] . The latter 
converges to the Oseen tensor for h — >■ oo. So, by continu- 
ously changing the height /i, one can control the range and 
strength of hydrodynamic interactions. They decay with 
l/r at large distances h from the wall, thus being of long 
range. In the vicinity of the wall, they decay with 
and justify the approximation of nearest-neighbor inter- 
actions used in Ref. [11]. In Fig. [8] in the inset we plot the 
relevant component Gfj of the Blake tensor as a function 
of h for, both, nearest-neighbor and next-nearest-neighbor 
rowers. The dashed lines indicate the respective values of 
the Oseen tensor that are reached in the limit h oo. 
On the other hand, for — >■ the relative strength of the 
nearest-neighbor component G22 compared to the next- 
nearest-neighbor component increases as the large graph 
in Fig. |8] demonstrates. So, we expect to observe with de- 
creasing h a transition from transient synchronization to 
the formation of metachronal waves. 
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Fig. 9: Average magnitude (a) and polar angle (b) of the 
order parameter Z for a rower chain situated at a height 
h above a bounding wall. The insets show the respective 
standard deviations. The global parameters of the open 
chain with = 200 rowers are /3 = tt/A, c ~ 0.28, a ~ 
—0.2, and e = —0.4. For h > 0.15, the total simulation 
time was T = 3750 and the time window for determining 
the averages was chosen sufficiently large. At — 2500, in 
order to minimize fluctuations of A and 3. For h < 0.15, 
the simulation time was extended to T = 6250 to ensure 
that a stationary metachronal state was reached and At = 
125. 



This is demonstrated in Figs.[9][a) and[9ljb), where we 
plot the time-averaged magnitude A and polar angle (p 
of the order parameter Z as a function of the height h. 
The insets depict the corresponding standard deviations. 
At hights h smaller than 0.35, A > 0.9 clearly indicates 
metachronism with an absolute phase difference between 
rowers smaller than 7r/4. It only varies slightly with h 
so that the wavelengths of the metachronal waves corre- 
spond to 9 to 10 rowers. Also the very small standard 
deviations are typical for metachronism. Then a relatively 
sharp transition occurs at h = 0.35 into a state of tran- 
sient synchronization. The amplitude A decreases by a 
factor of two and the negative ^ grows towards zero. In 
particular, the destruction of metachronism is indicated 
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Fig. 10: Color-coded phase differences of neighboring row- 
ers as a function of time in a chain of 400 rowers. The 
simulation starts at a height of /i = 0.20, which is then 
increased to ft- = 0.39 and h — 1.20 at times t — 500 
and t ~ 1000, respectively. At ft = 0.20 a metachonal 
wave develops and finally stretches over the entire chain. 
At ft = 0.39 the metachronism starts to break apart and 
transient synchronization with larger patches of correlated 
rowers are visible. At ft = 1.20 the patches become smaller. 



by a strong increase in the standard deviations showing 
that amplitude A and polar angle <P of the order parameter 
fluctuate strongly. In Fig. [TU] the transition is also visible. 
We show the color-coded phase differences between neigh- 
boring rowers as a function of time. At t = we start at a 
hight ft = 0.20 and then during the simulation switch to 
ft = 0.39 at t = 500 and to ft = 1.20 at t = 1000. Whereas 
the nearly uniform blue region at ft = 0.20 shows the 
existence of metachronal waves, they dissolve with time 
at ft = 0.39 and transient patches of same color indicate 
some synchronization. At ft = 1.20 the patches become 
smaller corresponding to a further decrease in A and the 
results are similar to our flndings in the unbounded fluid 
[compare green graphs in Figs.^a) and|3Jb)]. 

Our results show that metachronal waves in a chain 
of in-phase synchronizing rowers can only develop in the 
vicinity of a surface, where the range of hydrodynamic in- 
teractions is strongly reduced. Long-range hydrodynamic 
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interactions destabilize the waves since they lead to the 
formation of transient clusters synchronized in-phase. 



5 Conclusions 

In order to shed light on the essential features for obtain- 
ing metachronal waves, we implemented a chain of oscil- 
lators or rowers that are driven by an external non-linear 
force derived from a driving-force potential. Two rowers 
synchronize in antiphase when the curvature of the force 
potential is positive or in phase when the curvature is neg- 
ative. We concentrated on the latter case and introduced 
a suitable order parameter for identifying metachronal 
waves. Our investigations show that metachronal waves 
exist only when the range of hydrodynamic interactions is 
restricted either artificially or by the presence of a bound- 
ing surface as in any relevant biological system. The neg- 
ative curvature in the driving force potential leads to me- 
tachronal waves with wavelengths of 7-10 rower distances. 
Curiously, in Paramecium a wavelength including eight cil- 
ias is observevd [5] . In comparison, positive curvature that 
promotes antiphase synchronization in a rower pair only 
gives shorter wavelengths of around four rowers [B] . In or- 
der to obtain waves that travel either in one or the other 
direction along the chain, we break its left-right symmetry 
by tilting the rower segment with respect to the normal of 
the chain and distinguish between a fast power and slow 
recovery stroke. These measures are only effective when 
the rowers are sufficiently close to each other (c < 1). 

Metachronal waves and their defining properties, such 
as the wavelength and the propagational direction rela- 
tive to the direction of the power stroke, are an essen- 
tially two-dimensional phenomenon as Machemer pointed 
out [S]. Therefore, it will be important to extend current 
research towards models defined on two-dimensional lat- 
tices. This should give further insight into the fascinating 
properties of metachronal waves. Theoretical work in this 
direction was carried out on microfluidic rotors [551IM1I55] . 
However, amongst the great variety of beautiful patterns 
evolving, in particular, in Ref. |35| . metachronal waves as 
seen on ciliated surfaces cannot be observed. We have re- 
cently started to look at two-dimensional regular arrays 
of mainly stiff rods attached to a surface. We let them 
perform strokes where they move on a cone tilted with 
respect to the normal of the surface [H] . When we divide 
the beating cycle into a fast power and a slow recovery 
stroke, the latter acting when the rod is close to the sur- 
face, metachronal waves emerge. For a planar stroke see 
Ref. |45| . We will continue working on our system and 
also think about closed curved surfaces where topological 
constraints arise. 



A Mobilities for a fluid bounded by a plane 
wall 

When the fluid is bounded by an infinitely extended wall 
in the xy plane, the mobilities for point-like particles as- 
sume the following expressions. The self-mobility depends 



on the distance h from the bounding wall |43) : 



mm ('^m) 



Mo 




(7) 



and the Blake tensor G'^^(r) provides the cross mobihties 
= G'^^r.r, - r„) - GO^(r„, - f„) (8) 



where r„ ~ {xn,yn,h) is the location of the point force 
and r„ = (x„, j/„, — /i) is its image point piBT]. Using the 
Kronecker symbol 6\^^ the components of the so-called 
source doublet and Stokes doublet G^' read, respec- 
tively. 



(9) 
(10) 
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